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Atmospheric neutrino oscillations with DeepCore 

THE IceCube Collaboration^ 

^See special section in these proceedings 



Abstract: IceCube DeepCore can study atmospheric neutrino oscillations through a combination of its low energy reach, 
as low as about 10 GeV, and its unprecedented statistical sample, of about 150,000 triggered atmospheric muon neutrinos 
per year. With the diameter of the earth as a basehne, the muon neutrino disappearance minimum and tau neutrino 
appearance maximum are expected at about 25 GeV, which is considerably lower energy than typical IceCube neutrino 

events, but higher than the energies at which accelerator-based experiments have detected oscillations, disappearance 
and i^T appearance from neutrino oscillation can be measured in IceCube DeepCore. We present here the status of the 
newly developed low energy reconstruction algorithms, the expected experimental signatures, and the proposed approach 
for such neutrino oscillation measurements. 
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1 IceCube DeepCore 

IceCube is a cubic-kilometer neutrino observatory located 
at the geographic South Pole. IceCube construction be- 
gan in 2004 and was completed in December 2010. The 
complete detector consists of 86 strings deployed into the 
glacial ice, each of which consists of 60 Digital Optical 
Modules (DOMs) located between depths of 1450 m and 
2450 m. Seventy-eight strings are arranged on a hexagonal 
grid with an average 125 m horizontal spacing and 17 m 
vertical DOM spacing. The remaining 8 strings are more 
closely spaced in the center of the detector, with horizon- 
tal distances of 40 - 70 m and vertical DOM spacing of 
7 m. The 8 inner densely instrumented strings, optimized 
for low energies, together with the surrounding 12 IceCube 
standard strings, form the DeepCore inner detector (Fig. 1). 
These 8 inner strings have 10 DOMs located between 1750 
m and 1850 m in depth and 50 DOMs located between 
2100 m and 2450 m. The ice at depths between 1970 m 
and 2100 m, formed about 65,000 years ago [1], has a rel- 
atively short absorption length, and is known as the "dust 
layer". The 10 DOMs on each DeepCore string deployed 
above the dust layer help reject cosmic ray muons which 
are the major background to atmospheric neutrino studies. 
The ice below 2100 m has a scattering length about twice 



that of the ice in the upper part of the IceCube detector 
[2]. The lower 50 DeepCore DOMs are deployed in this 
very clear ice. DeepCore DOMs contain high quantum ef- 
ficiency photomultipliers (HQE PMTs [3]) which add ~ 
35% increase in efficiency compared to the standard Ice- 
Cube DOMs. With denser string and DOM spacing, clearer 
ice, as well as higher efficiency PMTs, DeepCore is opti- 
mized for low energy neutrino physics [4]. Fig. 2 shows 
the predicted and i/g effective areas at both trigger and 
onhne veto levels. Below 100 GeV the addition of Deep- 
Core increases the effective area of IceCube by more than 
an order of magnitude. 

2 Neutrino Oscillation Physics in DeepCore 

The IceCube DeepCore sub-array has opened a new win- 
dow on atmospheric neutrino oscillation physics with its 
low energy reach to about 10 GeV. The oscillation mea- 
surement is also made feasible by DeepCore' s location at 
the bottom center of IceCube, which allows the surround- 
ing IceCube strings to act as an active veto against cos- 
mic ray muons, the primary background to atmospheric 
neutrino measurements. A muon background rejection of 
8 X 10-3 for the overall IceCube trigger (~ 2000 Hz) is 



S.EULER etal. FUNDAMENTAL ATMOSPHERIC NEUTRINO MEASUREMENTS 



Overhead View 



OeepCore 



> 
(?) 



Hi 



• IceCube Strings 

HQE DeepCore Strings 
A DeepCore infili Strings 

(Mix of HOE and 

normai DOMs) 



DeepCore strings have 
10 DOMs with a 
DOiUI-to-DOM spacing 
of 10 meters 



' ' ' ' ' ' Dust Layer 



50 HQE DOMs with an 
DOM-to-DOM spacing 
of 7 meters 

21 Normai DOMs with a 
DOM-to-DOM spacing 
of 17 meters 



Figure 1: Overhead and side views of the IceCube Deep- 
Core detector The shaded hexagon in the overhead view 
shows the area covered by the DeepCore sub-array. On 
the side view, the hatched region shows the dust layer, and 
shaded boxes indicate the location of the DeepCore DOMs, 
10 in each string above the dust layer, and 50 in each string 
below the dust layer. 



achieved by applying a veto algorithm which rejects events 
with particle speed (defined as the speed of a particle trav- 
eling from the hit in the surrounding region to the center 
of gravity (COG) in DeepCore) between 0.25 m/ns and 0.4 
m/ns. The scheme of this veto algorithm is demonstrated 
in Fig. 3. With further expected improvements in the veto 
and event reconstruction algorithms, DeepCore expects to 
achieve a cosmic-ray muon rejection factor of 10^ or bet- 
ter [6]. Specific methods to investigate the oscillation phe- 
nomenon will be discussed in detail in the following sub- 
sections. 



2.1 Disappearance 

The earliest atmospheric neutrino oscillation evidence can 
be traced back to the zenith angle dependence of the dou- 
ble ratio measurement at few GeV energies in Super- 
Kamiokande [7]. IceCube DeepCore, with its approxi- 
mately 13 MT fiducial volume, is capable of making atmo- 
spheric neutrino oscillation measurements above 10 GeV, 
an energy region that has not been well explored by pre- 
vious atmospheric neutrino oscillation experiments. From 
Fig. 4, a significant deficit in the neutrino flux at 25 GeV 



is expected from the v^^ 



survival probability. In 



Fig. 5 this disappearance signature is shown for one year 
of simulated DeepCore data. The disappearance signal as- 
sumes that the path length is the diameter of the Earth, and 
therefore the ideal neutrino sample should contain only up- 
going neutrino-induced muons. An intrinsic difficulty for 
all experiments in identifying perfectly up-going neutrino- 
induced muon tracks is that the average opening angle, 
defined as the angle between the final state lepton direc- 
tion and the incoming neutrino direction, increases with 
decreasing energy. The uncertainty in the opening angle 
can be approximated as A* ~ 30° x .J {IGeV) / E^^. 
This effect will smear the oscillation signature in the neu- 
trino flux at lower energies. However, as shown in Fig. 5, 
DeepCore simulations indicate the potential to measure 
oscillations even if "up-going" tracks are defined to in- 
clude measured directions over the wide range of values, 
-1.0 < 008(6*) < -0.6 [4]. 

2.2 v-r Appearance 

The OPERA neutrino detector, located at the underground 
Gran Sasso Laboratory (LNGS), was designed for direct 
observation of — > appearance. OPERA announced a 
first tau lepton candidate from j/^ oscillation, and continues 
operation to achieve a statistically significant observation 
[8, 9]. DeepCore is currently acquiring data and will col- 
lect the world's largest inclusive sample of Vr- From Fig. 4, 
the region where flux reaches its minimum is the same 
region where Vt flux shows its corresponding maximum 
from z-'p — > Vt oscillation, that interact in DeepCore will 
produce an electromagnetic or hadronic shower, or "cas- 
cade". Events with short tracks which are beyond Deep- 
Core's ability to separate from cascades, are an irreducible 
background to cascade-like events. Therefore, cascade- 
like events include: 1) neutral current (NC) events from all 
three neutrino flavors (e, /i, r), 2) charged current (CC) 
events with short tracks (< ~' O (10) m), 3) Vr charged 
current events with r-leptons decaying into electrons or 
hadrons, 4) charged current events with r-leptons de- 
caying into muons whose track length is less than O (10) 
m. DeepCore should detect an excess of cascade-like 
events due to oscillation compared to the number of cas- 
cade events expected without oscillation. DeepCore may 
also be able to detect a distortion in the energy spectrum 
of cascade events due to v-r appearance. The simulated ex- 
cess of cascade-like events above ^ 25 GeV in DeepCore 
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Figure 2: Effective areas for muon neutrinos (left) and electron neutrinos (right). Triangles: IceCube standard strings 
only, trigger level. Squares: IceCube including DeepCore strings, trigger level. Circles: IceCube DeepCore after applying 
the online veto. 



is shown in Fig. 6. The deficit of oscillated cascade-like 
events compared to unoscillated below ^ 25 GeV is due to 
the rapid oscillation in the and ly-r survival probabilities 
shown in Fig. 4. Within the rapid oscillation regime, a large 
fraction of charged current (CC) events oscillate into 
events. The produced tau lepton from those i^t CC events 
always decay to at least one neutrino, reducing the visible 
energy and possibly resulting in an event below the detector 
energy threshold. However, without oscillation, these low 
energy i/^ CC events would be classified as cascade-like 
due to their short tracks. This significant deficit between 
the oscillated and unoscillated cascade-like event rates near 
the detector threshold may offer an opportunity to measure 
i/p disappearance in the cascade channel. The feasibility of 
measuring the excess of cascade-like events above 25 GeV, 
as well as systematic effects due to the ice properties, DOM 
efficiency and other factors, are under study. 

2.3 Reconstruction Algorithms 

Several reconstruction algorithms are being developed 
specifically for low-energy analysis in DeepCore. With 
shorter muon track lengths, low-energy events include a 
higher proportion of "starting" and fully contained muon 
events as opposed to through-going muons. The oscillation 
analysis also depends on separating track-like charged 
current events from cascade-like all-flavor neutral current 
and j/e and charged current events. Fully-contained- 
muon reconstruction algorithms calculate the length of the 
track from the reconstructed beginning and end points of 
the track. Fig. 7 shows the difference between recon- 
structed and true track length from simulated data recon- 
structed with one such algorithm. The starting muon event 
signature consists of a cascade associated with the charged- 
current muon neutrino interaction and a track associated 
with the resulting muon. Reconstructions of such events 




vertex time 
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^ particle speed 

Figure 3: Scheme of veto algorithm based on simulation. 
Upper part demonstrates a down-going muon hitting the 
DOMs with hit sequence varying from earliest in red to 
latest in blue. The big black circle exhibits the COG of 
the hits in DeepCore. Bottom part includes the vertex time 
and particle speed per hit. The "cut region" illustrates a 
cut based on the particle speeds which are consistent with 
down-going muons. 
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Figure 4: Oscillation probabilities for with sin^29i3 
0.1, Am§i = 2.5 X lO-^eV^,L = 12757 km [5]. 
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Figure 5: Simulated disappearance with one year Deep- 
Core data. Filled squares: Distribution of number of hit 
DOMs (NChannel) assuming no oscillation. Empty circles: 
Distribution of number of hit DOMs (NChannel) with os- 
cillation. 
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Figure 6: Simulated i^t appearance with one year Deep- 
Core data. Dotted lines: Distribution of true neutrino en- 
ergy assuming no oscillation. Solid lines: Distribution of 
true neutrino energy with oscillation. 
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Figure 7: Difference between reconstructed track length 
and Monte Carlo simulation true track length. 

therefore include the contributions of the cascade and the 
track. An algorithm under development for seeding these 
reconstructions and thereby separating tracks from pure 
cascades is based on a Fermat Surface [10]. 

3 Conclusions 

Simulations suggest that IceCube DeepCore may have the 
capability to measure atmospheric neutrino oscillations in 
an energy range which complements existing accelerator 
measurements. DeepCore construction is complete and the 
detector is already collecting data. New reconstruction al- 
gorithms suited to low energy measurements will enable 
IceCube to fully exploit the physics capabilities of Deep- 
Core. 
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Supernova detection with IceCube and beyond 

The IceCube collaboration^ 
^See special section in these proceedings 



Abstract: In its current configuration, IceCube is a formidable detector for supernovae. It can detect subtle features 
in the temporal development of MeV neutrinos from the core collapse of nearby massive stars. For a supernova at 
the galactic center, its sensitivity matches that of a background-free megaton-scale supernova search experiment and 
triggers on supernovae with about 200, 20, and 6 standard deviations at the galactic center (10 kpc), the galactic 
edge (30 kpc), and the Large Magellanic Cloud (50 kpc). Signal significances are reduced due to the noise floor and 
correlations between background hits. In this paper we discuss ways to improve the signal over background ratio with 
an improved data acquisition. We also discuss methods to track the average neutrino energy by multiple hit detec- 
tion from individual interacting neutrinos. The latter relies on the ability to reject coincident hits from atmospheric muons. 
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1 Introduction 

It was recognized early by [1] and [2] that neutrino tele- 
scopes offer the possibility to monitor our Galaxy for su- 
pernovae. IceCube is uniquely suited for this measurement 
due to its location and 1 km"^ size. The noise rates in Ice- 
Cubes photomultiplier tubes average around 540 Hz since 
they are surrounded by inert and cold ice with depth de- 
pendent temperatures ranging from —43 °C to —20 °C. At 
depths between (1450 — 2450) m they are partly shielded 
from cosmic rays. Cherenkov light induced by neutrino in- 
teractions will increase the count rate of all light sensors 
above their average value. Although this increase in indi- 
vidual light sensor is not statistically significant, the effect 
will be clearly seen once the rise is considered collectively 
over many sensors. 

The 5160 photomultipliers are installed in modules called 
digital optical modules (DOMs) and arranged in two con- 
figurations: IceCube, with 17 m (125 m) vertical spac- 
ing between DOMs (horizontal spacing between strings), 
and DeepCore with 7 m (72 m) spacings and equipped 
with high quantum efficiency DOMs, where eDocpCoro = 
1.35eicoCubo [3, 4]. With absorption lengths exceeding 
100 m, photons travel long distances in the ice such that 
each DOM effectively monitors several hundred cubic me- 
ters of ice. The inverse beta process z^e +p — > e+ +n domi- 
nates supernova neutrino interactions with 0(10 MeV) en- 



ergy in ice or water, leading to positron tracks of about 
0.55 cm • Ey/McV length for < 10 MeV. Considering 
the approximate dependence of the cross section, the 
light yield per neutrino roughly scales with E^. The detec- 
tion principle was demonstrated with the AMANDA exper- 
iment, IceCubes predecessor [5]. Since 2009, IceCube has 
been sending real-time datagrams to the Supernova Early 
Warning System (SNEWS) [6] when detecting supernova 
candidate events. 

Currently, the supernovae search algorithms are based on 
count rates of individual DOMs stored in 1 .67 ms time bins. 
We plan to introduce an improved data acquisition system 
that will allow to store all IceCube hits for supernova candi- 
dates. We discuss below some of the improvements that we 
expect to be achieved using this additional information. In 
addition, the collaboration is discussing future extensions 
of the detector that would also improve the supernova de- 
tection capacity. 

2 Current performance 

IceCube is the most precise detector for analyzing the neu- 
trino lightcurve of close supernovae. A paper, discussing 
the detector and physics performance, is close to being pub- 
lished. Figure 1 shows the expected significance for the 
detection of a supernovae as function of distance (left) and 
presents the expected rate distribution for the Lawrence- 
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Figure 1 : Left: Expected significance versus distance assuming the Lawrence-Livermore model [9] for three oscillation 
scenarios. The significances are increased by neutrino oscillations in the star by typically 40% in case of an inverted 
hierarchy. The Magellanic Clouds as well as the center and the edge of the Milky Way and various trigger thresholds 
are marked. For the Milky Way, the supernova progenitor distribution follows the prediction from [7], for the Magellanic 
Clouds it is assumed to be uniform. Right: Expected rate distribution at lOkpc supernova distance assuming normal and 
inverse hierarchies. 
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Figure 2: Left: Measured probability density distribution of time differences between pulses for noise (bold line) and the 
expectation for a Poissonian process fitted in the range 15 ms < AT < 50 ms (thin line). The excess is due to bursts 
of correlated hits, as indicated by the 50 ms long snapshot of hit times shown in the inset. Right: Measm^ed width of 
the significance distribution as function of time during IceCube construction with 40 (left), 59 (middle) and 79 (right) 
deployed strings. The inset shows the significance distribution before (wide distribution) and after (narrow distribution) 
suppression of hits due to atmospheric muons (79 strings). The current trigger threshold for SNEWS alarms is indicated 
by a dashed line. 



Livermore model [9] (right). The rate distribution demon- 
strates the excellent resolution of details in the neutrino 
light curve. This includes the possibility to distinguish the 
neutrino hierarchy, provided that the astrophysical model 
is well known and sin^ 813 > 1°[7]. The present noise 
floor is indicated in Fig. 1 (right) which leads to a fast de- 
terioration of the signal significance particularly at larger 
distances. In addition, the expected signal significance in 
IceCube is somewhat reduced due to two types of correla- 
tions between pulses that introduce supra-Poissonian fluc- 
tuations. The first correlation involves a single photomul- 
tiplier tube. It comes about because a radioactive decay 
in the pressure sphere can produce a burst of photons last- 
ing several /xs. The second correlation arises from the cos- 
mic ray muon background; a single cosmic ray shower can 



produce a bundle of muons which is seen by hundreds of 
DOMs. The observed time difference between noise hits 
deviates from an exponential distribution expected for a 
Poissonian process (see Fig. 2, left). The inset shows a 
hit sequence from a single DOM, clearly indicating the 
bursting behavior. A significant fraction of these bursts 
can be rejected by an artificial non-paralyzing deadtime, 
currently adjusted to t = 250 /is, which decreases the av- 
erage optical module noise rate to 285 Hz, while keeping 
« 87%/(l + rsN ■ ''') of supernova induced hits with rate 

Due to remaining correlated pulses from radioactive decays 
and atmospheric muons, the measured sample standard de- 
viation in data taken with 79 strings is « 1.3 and sa 1.7 
times larger than the Poissonian expectation for 2 ms and 
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IceCube single hits, eff. V = 390 m^ 
[ceCube NN hits, eff. V = 0.548 m^ 
DeepCore single hits, eff. V = 527 m^ 
DeepCoreNN hits, eff. V = 2.417 m^ 
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Figure 3: Effective volume per optical module and positron 
energy versus positron energy. Due to the chosen generation 
volume, the single hits effective volumes (circles) are about a 
factor of two too small. This is illustrated by the calculations 
taken from [11] (lines). The nearest neighbor coincidence 
modes (squares) are shown here without cut on the coinci- 
dence time. The effective volumes given in the legend are 
calculate according to the "Sf model from [12] with an in- 
tegration time of 4 s. 
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Figure 4: Left: Significance as a function of distance d and the allowed time window [25 ns — t, 25 ns + t] between hits in 
two neighboring DOMs 7 m apart (DeepCore). The time difference distribution is shown in the inset in solid for DeepCore 
and dashed for IceCube. The maximum distance of about 10 kpc is reached for i = 25 ns for a significance threshold of 
6.7 (white line). Right: Energy resolution AE as a function of distance d and average energy of the emission spectrum 
(E). The white line show several AE/E contours. 



500 ms time bins, respectively. For the analysis, we de- 
fine the significance f = Aji/a/^^ , where A/^t is the most 
likely collective rate deviation of all optical module noise 
rates from their running average, a^^ is the correspond- 
ing deviation calculated from the data, thus accounting for 
non-Poissonian behavior in the dark rates. The signifi- 
cance should be centered at zero with unit width. Fig. 2 
(right plot) shows that this is not the case; instead the width 
changes with the season and increases with the size of the 
instrumented detector This behavior is linked to the sea- 
sonal change of the muon flux. While atmospheric muons 
contribute to the count rates of individual DOMs by only 
3%, these hits are correlated across the detector, thus broad- 
ening the significance distribution and giving rise to a sea- 
sonal dependence of the trigger rate. At present, it is possi- 
ble to subtract roughly half of the hits introduced by atmo- 
spheric muons from the total noise rate offline, as the num- 
ber of coincident hits in neighboring DOMs is recorded for 
all triggered events. The width of the significance distribu- 
tion then decreases to about 1.06, close to the expectation 
(see inset of Fig. 2 (right plot)). 

A data acquisition that records all hits in case of a super- 
nova trigger will permit further improvements. The time 
resolution on the onset of the burst will no longer be re- 
stricted by the 1.67 ms time bins in which the rates are 



recorded, hits associated to triggered atmospheric muons 
can be fully rejected, and more sophisticated methods to 
minimize correlated pulse bursts, e.g. by eliminating the 
initial hits of the bursts while keeping photomultiplier re- 
lated afterpulses, can be applied. However, as the signifi- 
cance improves only with l/-\/A^background, a much more 
drastic reduction of background is required in order to im- 
prove the detection significance at the edge of our galaxy, 
to track the average neutrino energy and maybe even pro- 
vide some directional capability. This can only be achieved 
by detecting more than one Cherenkov photon from an in- 
teraction and applying a coincidence condition. 

3 New opportunities from coincidence rates 

The study in this section is motivated by an analytical 
framework that explores the potential of coincident hit 
modes [11]. Here, we investigate the "nearest neighbor co- 
incidence hits" mode, with a hybrid GEANT-4/toy Monte 
Carlo simulation. We chose this mode from other possible 
multi-hit modes such as multiple hits in one DOM or co- 
incident hits between any DOMs, because it has the best 
noise suppression potential by requiring a very short time 
window around the two coincident hits. 
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Figure 3 shows the positron effective volume per module 
and positron energy for the two detector configurations and 
two detection modes: single hits and nearest neighbor co- 
incidence hits. The effective volumes given in the legend 
are calculated according to the "Sf model from [12] with 
an integration time of 4 s. Both modes have a very different 
energy dependence, which makes the ratio A'coinc /^single 
an observable of the average energy {E) of the emission 
spectrum. 

The inset in Fig. 4 shows the smallest time difference be- 
tween hits in neighboring DOMs for both detector configu- 
rations. A cut on this time distribution was found by calcu- 
lating the detection significance as a function of the super- 
nova distance and the time window, as is shown in Fig. 4 
for DeepCore. Applying a time cut of ±25 ns around the 
most probable time difference of To = 25 ns, a maximum 
distance of 10 kpc can be reached with the DeepCore DOM 
separation. A similar cut yields a smaller reach for IceCube 
due to the larger DOM separation. 

To estimate the energy resolution of the energy observable 
-^coincZ-^singie, the ratio was calculated using spectra ac- 
cording to the "Sf model with average energies between 
2 - 30MeV. In Fig. 4 (right plot), the deviation Ai? is 
shown as a function of the supernova distance and the av- 
erage {E) . For distance smaller than the trigger threshold 
(w 10 kpc), a resolution of around 5% can be achieved for 
spectra with an average energy of 10 — 15 MeV. 

The energy resolution depends on the expected noise level 
for the chosen selection mode. Above, the Poissonian noise 
levels were scaled up by 1.3, as mentioned in Sec. 2. It 
is possible that this underestimates the average hit prob- 
ability for DOMs that were close to atmospheric muons. 
Preliminary studies show that light from muons can be sup- 
pressed by considering only DOMs that were around 300 m 
or further away from the reconstructed track position. In 
the worst case, when a track traverses the whole detector 
volume vertically, this cut reduces the usable volume by 
« 30%. Alternatives to such a cut are being investigated. 

4 Possible extensions of IceCube/DeepCore 

Discussions on an extension of IceCube/DeepCore have 
started, which would also improve the supernova detection 
capability. We used a GEANT-4 simulation to estimate the 
capabilities of a hypothetic 18 string detector with IceCube 
DOMs spaced apart a few meters. Fig. 5 shows the effec- 
tive volumes of single and multiple hits as function of dis- 
tance between the DOMs. The sizable increase in the active 
volume of coincidence hits would strongly improve the sig- 
nal over noise ratio and lead e.g. to a substantial improve- 
ment of a supernova detection at the Magellanic cloud. 

5 Conclusion 

As a supernova detector, IceCube already offers an un- 
matched ability to establish subtle features in the tempo- 
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Figure 5; Effective volumes per DOM for multiple hits 
as function of horizontal and vertical distance between 
DOMs. The same neutrino spectrum was used as in Fig. 3 

ral development of the neutrino flux by tracking the overall 
count rates of its DOMs. The correlated noise background 
remains one of the big challenges. The new data acquisi- 
tion will permit the recording of all hits in case of a trigger 
signal, thus greatly improving the rejection of correlated 
noise sources such as atmospheric muons. It will also allow 
to study the influence of the correlated noise on multiple hit 
detection modes. 

One such mode, nearest neighbor coincident hits, was in- 
troduced in this paper and shows great potential to extent 
the existing capabilities of IceCube by measuring the aver- 
age neutrino energy. This mode also has the potential to be 
sensitive to the neutrino direction. Other multi-hit modes 
will be studied next. All these modes would greatly bene- 
fit from a very dense sub-array with inter-DOM and inter- 
string distance of a few meters only, extending their reach 
to several tens of kilo parsec and possibly beyond. 
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Abstract: The IceCube observatory, 1 km in size, is now complete with 86 strings deployed in the antarctic ice. 
IceCube detects the Cherenkov radiation emitted by charged particles passing through or created in the ice. To realize 
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IceCube data with simulation based on the new model is shown. 
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1 Introduction 

The properties of photon propagation in a transparent 
medium can be described in terms of the average distance 
between successive scatters and the average distance to ab- 
sorption (local scattering and absorption lengths), as well 
as the angular distribution of the new direction of a photon 
relative to old at a given scattering point. These details are 
used in both the simulation and reconstruction of IceCube 
data, thus they must be known to the best possible accuracy. 
This work presents a new, direct fit approach to determine 
these ice properties, which is different from the method de- 
scribed in [4]. A global fit is performed to a set of data 
with in-situ light sources (see Figure 1) covering all depths 
of the detector, resulting in a single set of scattering and ab- 
sorption parameters of ice, which describes these data best. 
Figure 2 shows examples of experimental data used for this 
analysis. 

2 Flasher dataset 

In 2008, IceCube consisted of 40 strings as shown in Figure 
3, each equipped with 60 equally spaced optical sensors, 
or digital optical modules (DOMs). Each of the DOMs 
consists of a 10" diameter photomultiplier tube (PMT) [2] 
and several electronics boards enclosed in a glass container 
[3]. One of the boards is the "flasher board", which has 
6 horizontal and 6 tilted LEDs, each capable of emitting 
~ 7.5 • 10^ photons at ^ 405 ± 5 nm in a 62 ns-wide pulse. 




Figure 1: Simplified schematics of the experimental setup: 
the flashing sensor on the left emits photons, which propa- 
gate through ice and are detected by a receiving sensor on 
the right. 

The PMT output signal is digitized into "waveforms" using 
the faster, ATWD, and slower, fADC, sampling chips [1]. 
The ATWD is configured to collect 128 samples with 3.3 ns 
sampling rate, and the fADC records 256 samples with 
25 ns sampling rate. The DOMs transmit time-stamped 
digitized PMT signal waveforms to computers at the sur- 
face. 

In a series of several special-purpose runs, IceCube took 
data with each of 60 DOMs on string 63 flashing in a se- 
quence. For each of the flashing DOMs at least 250 flasher 
events were collected and used in this analysis. All 6 hori- 
zontal LEDs were used simultaneously at maximum bright- 
ness and pulse width settings, creating a pattern of light 
around string 63 that is approximately azimuthally sym- 
metric. 
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Figure 2: Example photon arrival time distributions at a 
sensor on one of the nearest strings (122 m away), and on 
one of the next-to-nearest strings (217 m away). Dashed 
lines show data and solid lines show simulation based on 
the model of this work (with best fit parameters). The goal 
of this work is to find the best-fit ice parameters, which 
describe these distributions as observed in data simultane- 
ously for all pairs of emitters and receivers. 
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Figure 3: IceCube 40-string configuration as operated in 
2008. String 63 (of DOMs that were used as flashers) is 
shown in black. IceCube parts installed in the following 
years (2009, 2010 as shown in the figure) lie in regions 
indicated approximately with dashed lines. 



The pulses corresponding to the arriving photons were ex- 
tracted from the digitized waveforms and binned in 25 ns 
bins, from to 5000 ns from the start of the flasher pulse 
(extracted from the special-purpose ATWD channel of the 
flashing DOM). To reduce the contribution from saturated 
DOMs (most of which were on string 63 near the flashing 
DOM) [2] the photon data collected on string 63 was not 
used in the fit. 



3 Ice parametrization 

The ice is described by a table of parameters 5e(405), 
a(405), related to scattering and absorption at a wavelength 
of 405 nm at different depths. The width of the vertical ice 
layers (10 m) was chosen to be as small as possible while 
maintaining at least one receiving DOM in each layer. Co- 
incidentally it is the same as the value chosen in [4]. 

The geometrical scattering coefficient b determines the av- 
erage distance between successive scatters (as 1/6). It is of- 
ten more convenient to quote the effective scattering coeffi- 
cient, be = b ■ {1 — (cos 9)), where 6 is the deflection angle 
at each scatter, () denote the expectation value. The absorp- 
tion coefficient a determines the average distance traveled 
by photon before it is absorbed (as 1 / a). 

4 Simulation 

The detector response to flashing each of the 60 DOMs on 
string 63 needs to be simulated very quickly, so that simu- 
lations based on many different sets of coefficients fee (405) 
and a(405) could be compared to the data. 

A program called PPC (photon propagation code [7]), was 
written for this purpose. It propagates photons through 
ice described by a selected set of parameters fee (405) and 
a(4G5) until they hit a DOM or get absorbed. No special 
weighting scheme was employed, except that the DOMs 
were scaled up in size (a factor 5 to 16, depending on the 
required timing precision), and the number of emitted pho- 
tons was scaled down by a corresponding factor (5^^ — 16^). 

The probability distribution f{9) of the photon scattering 
angle 9 is modeled by a linear combination of two functions 
commonly used to approximate scattering on impurities: 

f{9) - (1 - /sl) • HG + fsL ■ SL, 

where HG is the Henyey-Greenstein function [4]: 



p(cos 6) 



1 



g = {cos t 



2[l + g2_2c,.cos6']3/2' 
and SL is the simplified Liu scattering function [8]: 

2.9 



p{cos 9) ^ {1 + cos 9)°' , with a = 



1-5 



/sL determines the relative fraction of the two scattering 
functions and it determines the overall shape. Figure 4 
compares these two functions with the prediction of the 
Mie theory with dust concentrations and radii distributions 
taken as described in [4]. The distributions of photon ar- 
rival time are substantially affected by the "shape" param- 
eter /sL (as shown in Figure 5). /sl is also a global free 
parameter in the fitting procedure. 

The value of g — 0.9 was used in this work (cf. g — 0.8 
in [4]). Higher values (as high as ^ 0.94 [4, 6]) are pre- 
dicted by the Mie scattering theory, however, these result 
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Figure 4: Comparison of the Mie scattering profiles cal- 
culated at several depths of the South Pole ice with the 
Henyey-Greenstein (HG) [4] and simplified Liu (SL) [8] 
scattering functions, all with the same g = 0.943. 
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Figure 5: Photon arriving time distributions at a DOM 
125 m away from the flasher, simulated for several values 
of g = {cos 9) and fsi- The difference in peak position 
simulated with g — 0.8 and g — 0.9 is of the same or- 
der (^10 ns) as that between sets simulated with different 
values of the shape parameter fsi- 



in slower simulation, while yielding values of the effec- 
tive scattering be and absorption a coefficients that change 
by less than 3% as determined in [4], which could also be 
concluded from Figure 5. 
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Figure 6: Values of 6e(405) and a(405) vs. depth for subse- 
quent steps of the minimizer. The "converged" black curve 
shows fitted values after the last of the 20 steps of the min- 
imizer 

and absorption in 10 m layers at detector depths of 1450 to 
2450 m). 

The photon counts d{ti) and s{ti) observed in time bins in 
rid data and Ug simulated flasher events are compared to 
each other using a likelihood function 



\/2ttc 



■ exp 



rf! 

-log Ms 



• R. 



The product is over all emitters and all time bins of re- 
ceivers. ^ld{ti) and iis{ti) are the expected values of pho- 
ton counts per event in data and simulation, and are deter- 
mined by maximizing C with respect to these. The first two 
terms in the product are the Poisson probabilities, and the 
third term describes the systematic uncertainties inherent 
in the simulation. The last term R represents regularization 
constraints of the solution values with depth and with each 
other. 

Starting with the homogeneous ice described with 
5e(405) 0.042 m-i and a(405) = 8.0 kvcT^ (average 
of [4] at detector depths) the maximum of C is found in ~ 
20 steps. At each iteration step the values of 6e(405) and 
a(405) are varied in consecutive ice layers, one layer at a 
time. Five flashing DOMs closest to the layer, which prop- 
erties are varied, are used to estimate the variation of the C. 
Figure 6 shows ice properties after each of 20 steps of the 
minimizer The general agreement of the model and data is 
good as shown in Figure 2. 



5 Fitting the flasher data 

Data from all pairs of emitter-receiver DOMs (located in 
the same or different ice layers, altogether ^ 38700 pairs) 
contributed to the fit of ~ 200 ice parameters (scattering 



6 Dust logger data 

Several dust loggers [5] were used during the deployment 
of seven of the IceCube strings to result in a survey of the 
structure of ice dust layers with extreme detail (with the 
effective resolution of ^ 2 millimeters). These were then 
matched up across the detector to result in a tilt map of the 
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Figure 7: Comparison of the average dust log with the 
effective scattering coefficient fee(405) measured with the 
flasher data. 

South Pole ice, as well as a high-detail average dust log, a 
record of a quantity proportional to the dust concentration 
vs. depth. Additionally, the EDML (East Dronning Maud 
Land, see [5]) ice core data was used to extend the dust 
record to below the lowest dust-logger-acquired point. 

The correlation between the effective scattering coefficient 
measured with the IceCube flasher data and the average 
dust log (scaled to the location of string 63) is excellent, as 
shown in Figure 7. Within the depth range 1450 m - 2450 
m instrumented with DOMs all major features match, have 
the right rise and falloff behavior, and are of the same mag- 
nitude. Some minor features are washed out in the flasher 
measurement. 

Having established the correlation with the average dust 
log, the EDML-extended version of the log was used to 
build an initial approximation to the fitting algorithm de- 
scribed in the previous section. This resulted in a solution 
that is determined by the scaled values of the extended log 
(instead of by the somewhat arbitrary values of the initial 
homogeneous ice approximation) in the regions where the 
flasher fitting method has no resolving power, i.e., above 
and below the detector. 



7 Results 

The effective scattering and absorption parameters of ice 
measured in this work are shown in Figure 8 with the 
±10% gray band corresponding to ±la uncertainty at most 
depths. The uncertainty grows beyond the shown band at 
depths above and below the detector. The value of the scat- 
tering function parameter fsi = 0.45 was also determined. 
Figure 8 also shows the AHA (Additionally Heterogeneous 
Absorption) model, which is based on the ice description of 
[4] extrapolated to cover the range of depths of IceCube and 
updated with a procedure enhancing the depth structure of 
the ice layers. 
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Figure 8: Values of 6e(405) and a(405) vs. depth for con- 
verged solution shown with solid lines. The updated model 
of [4] (AHA) is shown with dashed lines. The uncertainties 
of the AHA model at the AMANDA depths of 1730 ± 225 
m are ~ 5% in and ^ 14% in a. The scale and numbers 
to the right of each plot indicate the corresponding effective 
scattering 1/be and absorption 1/a lengths in meters. 
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